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ABSTRACT 


Observed ocean acoustic wavefronts show surprising stability in long-range 
acoustic transmission experiments. This suggests that ocean scattering processes tend to 
redistribute acoustic energy dominantly along the wavefront rather than across it. The 
purpose of this thesis is to elucidate the physical mechanism for this type of scattering by 
presenting a ray-based physical model. Analytic formulae are presented that predict 
wavefront distortions caused by ocean internal waves and other processes. Applications 
of this study include wavefront healing near underwater obstacles, out-of-plane scattering 


and the vertical redistribution of energy of off-axis sources. 
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I. INTRODUCTION 


Observed ocean acoustic wavefronts show surprising stability in long-range 
acoustic transmission experiments. This suggests that ocean scattering processes tend to 
redistribute acoustic energy dominantly along the wavefront rather than across it. The 


physical mechanism for this type of scattering has not yet been elucidated. 


There are several practical implications for this type of scattering. First, 
scattering along the wavefront leads to the important result that shadow zones behind 
ocean obstacles, like seamounts, will be filled in more rapidly than by simple diffraction. 
Second, this type of scattering can lead to significant transfer of acoustical energy from 
high angles to low angles. Hence, an off-axis source, like those responsible for the ocean 
ambient sound fields, will have their energy redistributed vertically, thus altering the 
vertical directionality of the ambient sound. Thirdly, scattering along the front can 
contribute to the vertical extension of acoustical caustics zones that abut shadow zones. 
Observations reported by Dushaw et al. (1999) of the Acoustic Thermometry of Ocean 
Climate group on Navy SOSUS stations showed that wavefront data obtained from these 
stations had experienced dramatic vertical extensions of wavefront caustics several 
hundreds of meters above the receivers. The resultant scattering caused interference 


patterns along the front and vertically extended the front into shadow zones (Figure 1). 


Finally, since this scattering mechanism leads to acoustical energy redistribution 
along the front, the scattering processes can be considered to be coherent in some regards. 
Thus, the scattered energy can be used effectively for signal and array processing needs, 


but the statistics from this type of scattering has not yet been examined in detail. 


With this background information in mind, we consider a pulse emitted from a 
point source in the ocean sound channel. In this case, an outward propagating wavefront 
is formed. After traveling approximately 1,000 km, the wavefront takes on an accordion 
shape (Figure 2). The wavefront is composed of individual segments that are joined by 
cusps and each point on the wavefront is identified with a ray with a specific launch angle 


of the source. The wavefront is thus defined by the end points of all the rays at a specific 


time after emission of the pulse. After adding small-scale heterogeneities due to internal 
waves, each ray will be displaced some distance from its unperturbed end point. The 
deviation of this ray end point can be either along the wavefront or perpendicular to the 
wavefront. Numerical simulations in this thesis and in other work will show that the 


deviation is much larger along the wavefront. 


The primary purpose of this thesis is to study along the front scattering using a 
hierarchy of models from a simple 1-D, discrete model to a fully random ocean model. 
The main acoustical observable is the geometry of the acoustical wavefront, and Ray 
methods are used to quantify the re-distribution of ray energy along and across the front. 
Of principal importance is the ability to predict the ocean and acoustic variables that 
could most affect this type of scattering. For example, many of our simple reduced 
physics models predict that variance of along wavefront ray deflections scale like the 
cube of the propagation range, while the variance of across the front ray deflections scale 
linearly with range. In addition, scattering in the wavefront direction is controlled by the 
vertical gradients of sound speed while the scattering in the ray direction depends on the 


sound speed perturbation itself. 


U. METHODOLOGY 


Ray-based acoustic propagation models with various ocean sound speed 
fluctuation models were used to estimate scattering along and across the wavefront. The 
simplest model considered was a discrete 1-D model with small angle forward scattering. 
Here Snell’s law was linearized and used to analytically predict the trajectory of ray paths 
through a series of random sound speed layers. The validity of the linear approximations 
was examined to determine the relative contributions of along vs. across wavefront 
scattering as a function of the sound speed variance, the number of layers, the width of 


the layers, and the wavefront angle of incidence. 


The next level of complexity considered a continuous 2-D random field of sound 
speed perturbations obeying the (GM) internal wave spectrum. Various analytic 
approximations to the scattering were checked through the use of Monte-Carlo numerical 
ray simulations that examined the scattering physics. The simplest case modeled was ray 
propagation through an internal wave field that was anisotropic in depth and the 
horizontal direction, but homogeneous in depth. A comparison of ray simulation 
statistics to some simple analytic expressions was easily obtained because of the simple 


geometry. 


The last and most complicated case involved a true ocean environment, with a 
sound channel, and an anisotropic and inhomogeneous internal wave field. Again, 
Monte-Carlo ray simulations were performed to test more complicated analytic 


expressions for scattering along and across the front. 


A. SOUND SPEED ENVIRONMENT 


In this numerical experiment, the following two dimensional sound speed field is 


considered: 


e(x,z)=C(z)+8c(x,z), dc c(z) (2.1) 














where dc= (= Cw and c(z) represent the sound speed perturbation and mean sound 
Zz 


Pp 


speed profile respectively. € is the internal wave displacement field and (= is the 
Z 
p 


potential sound speed gradient. The mean sound speed profile can then be modeled in 
one of two ways. Either, it is represented by the constant co=1500 ms"! or by the more 


complex Munk canonical form (Munk, 1974). 


G(z) =co| I+e(e%)"" +2(z-m)/B-1) (2.2) 


Specifically, the axial sound speed (co) is equal to 1500 ms", the total water depth 
D=5000 m, B=1000 m, the sound channel axial depth (z,) was set to 1000 m, and 
¢ =0.005515. 


First, the potential sound-speed gradient profile is expressed as 


h(nt)= te), [=] 2.3) 


where Co is 1500 ms", N is the buoyancy frequency and a is a dimensionless constant 
yancy Ireq y 


related to salinity and temperature, approximately equal to 24.5. 


From Flatte et al (1979), a fractional sound speed variance profile was imposed in order 
to understand the internal wave induced sound-speed fluctuations: 


(H° (z)) =(u,) ne) (2.4) 


3 
Nit 


where (u,) =6.26 x 10°, and Nret = 3 cph. However, near the surface and bottom, the 


profile of (Ww? (z)) will be modified in order to accommodate the internal wave vertical 


mode boundary condition of zero displacement at those areas. Equation 2.4 is important 
because it reveals the inhomogeneities of internal wave induced sound-speed fluctuations 
with respect to depth. As a result of a larger buoyancy frequency in the upper ocean as 


opposed to a deeper one, these sound-speed fluctuations are larger in the upper regime. 


4 


The task now falls to generating the random realizations of sound speed 
perturbations that are induced by internal waves, where, internal waves are oscillatory 
geophysical wave motions resulting from the Coriolis and buoyancy forces. The change 
in sound-speed due to internal waves is on the order of 1 ms”, with time scales lying 
within the regime of the Coriolis and buoyancy periods. The length scales of internal 
waves are of order 100 m in the vertical and 10 km in the horizontal. Further, sound- 
speed fluctuations caused by internal waves are the dominant source of variability of 


acoustic fields in the ocean. 


The Garrett-Munk internal wave spectrum is utilized using the Colosi and Brown 
method (1998). Internal waves can be represented as horizontally propagating vertical 


modes. Consequently, the internal wave displacement field ¢ can be considered as the 
sum over all possible internal waves with vertical mode number j and wave numbers kx 
and ky. 


Jmax 


Cw (%Y.Z) = Re i dk, })g (jk, )A, (ze (2.5) 


oe jel 


where Aj is the a mode function, k, is the horizontal wave number in the propagation 


plane. The amplitudes of the internal waves are complex random variables with statistics 
described by the GM spectrum. The buoyancy frequency profile is modeled using a 


canonical profile as shown below. 


N(z) =Nse”” (2.6) 


Here, B is still equal to 1,000 m and the surface buoyancy frequency (Ns) is set to 5 cph. 
However, when the background sound speed is a constant, N.=Icph is employed. The 
sum of the contributions from internal wave mode numbers up to jmax = 200 was chosen 
and internal waves with horizontal wave scales from 0.5 to 409.6 km were utilized in the 
generation of the random realizations of internal wave induced sound speed 


perturbations. 


A random realization of the spectrum of 1 is achieved by expressing it in terms 


of a mode number and horizontal wave numbers kx and ky. The final result is 
5 


gN el 2) Ty ee (ki+k;)? +k, 
3 1 


TN MPT EAT 2 iayey | (cae)! ok 





G,(j,k,) =H (2.7) 


fj . a ; Bh ace : : 
where k, = “71 | Internal waves are isotropic in the horizontal direction, but anisotropic 
0 


in the vertical. 


B. RAY THEORY 


The evolution of the wavefront is well described by the two-dimensional ray 
equations (Brown et al, 2003, and Beron-Vera et al 2003). They are given by 


dx >.2\/2 dz 4, dp lé 
—=c(l-c = =< = 2.8 
dt ae a at ce ee 


where p=sin@/c is the ray slowness in the z-direction. The ray equations were integrated 
forward in time to Tinax=1000 km/co = 666.667(s), with the intermediate times being 


saved in increments of 0.667 seconds. These equations were solved numerically using 


MATLAB. 


C; NO WAVEGUIDE CASE 


For the case of constant background sound speed, simulations were carried out 
with a source depth of 1,000 m and with an initial ray angle of zero degrees. A total of 
12,000 realizations of ray propagation were computed in order to accurately calculate the 
necessary statistics. In this case, a periodic boundary condition was implemented in 
order to compensate for rapid scattering to the ocean boundaries. Figure 3 shows the ray 
endpoints from all the realizations at the following nominal ranges: 340, 680, and 1,000 
km. A dramatic depth deviation (wavefront direction) of the end point relative to the 
range deviation (ray direction) is clearly evident. This figure dramatically shows the 
dominance of scattering along the front. The rays are also seen to be scattered into an arc 
instead of a straight line. The arc is formed because the rays with the largest depth 
deviation could not get down range as far as the rays that stayed closer to the original 


depth. In Figure 3, the green curve represents the parabolic fit to ray end point 
6 


distributions and the red curve shows the location of the unperturbed point source front 
(spherical wave). The parabolic fit is used as a mean wavefront and thus this function is 
subtracted off the ray endpoints to yield statistics in the wavefront and ray directions. It 
is clear that at 1000 km the wavefront direction deviations are on the order of tens of 
kilometers as opposed to the order of 10 meters of the ray direction deviations. Statistical 


results will be given in later sections. 


D. WAVEGUIDE CASE 


Implementing a source depth of 1,000 m and an initial ray launch angle 


of 6:(0) =+7°, where the angle @: is with respect to the x-axis, the case of the Munk 


background profile with internal wave induced sound speed variations was simulated. A 
total of 12,000 realizations of ray propagation were modeled. Figure 4 shows the ray 
endpoints from all the realizations at nominal ranges of 340, 680, and 1000 km. The 
location of the unperturbed +7° ray is denoted by the green dots. Further, the 
unperturbed front composed of initial angles between +0° and +14° is denoted with a red 
line. At a range of 1000 km, deviations of the ray end points in the ray direction are of 
the order 10 m. Again, in the wavefront direction, the order is in kilometers, 
demonstrating the dominance of along wavefront deviations. Also note that the 
scattering in the wavefront direction is seen to be asymmetrical in the sense that the 
wavefront is extended more towards the longer and not the shorter ranges. Statistics will 


be discussed in later sections. 
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Hl. 1-D DISCRETE CASE 


The oceanographic wave propagation physics leading to distinct anisotropy in 
wavefront distortions is geometrically complex. To alleviate this complexity, and better 
understand the precise geometry involved with this problem, a plane 1-D model is first 
considered. Although, the model is not realistic, its simplicity allows for the analytical 
solutions of the trajectories of individual rays which are necessary in elucidating the 


aforementioned geometry. 


A model consisting of a layered random sound speed medium with small angle 
forward scattering was considered (see Figure 5). There is an incident wave with angle 
Oo (with respect to the horizontal), in a region of sound speed co. The incident wave 


refracts via Snell’s Law with angular deviation 56: into region 1| of a slab with specified 











width L, where Colosi and Flatte (2007) defined the sound speed asc,=c,+6c,,6c, 1] Cy. 





The scattered ray then moves a distance x; in the horizontal, taking time t; to cross the 
slab in this region. The process continues in this manner for a specified number of N 


slabs. 


Colosi and Flatte (2007) used Snell’s Law, cos? 
c 





= constant, to express a small 


angular deviation when there is a small sound speed perturbationédc,;. Their results are 


given by 





§6;z—3 (3.1) 


where p,;=6c,/c, and@,=@,+060,. After the scattered ray crosses each layer, an x- 


deflection and a travel time result. Their resulting equations are exact relations. 


(2 ee 3.2) 
tan(@0 + 6)) cj sin(@0+ 6)) 





After expansion in small angles, Flatte and Colosi (2007) arrived at the following 


approximation: 


cos 200 
,oti= = ti- -to= top 
sin Be sin’ Go 








OXj = Xj-Xo = Xo (3.3) 

Here, the unperturbed deflections of the ray path are represented by the values 
xo=L/tan@o andt,=L/(c,sin@,). According to the results (Equation 3.3), both the x- 
deflection 5x; and the sound speed perturbation pj will always have the same sign. That 
is to say, should pj be positive then 560; must be negative, which results in the ray 
having a larger glancing angle and thus a largerdx,;. The opposite also holds true. 
However, the sign of the time shift is dependent upon both the angle of incidence and the 
sound speed perturbation. For the case of interest @0 < 45°, the time shift dt; always 
holds the same sign as the x-deflection. With pj > 0, the path length is longer and the 


time shift positive. Thus, the sound speed increase is unable to compensate for the 
increase in path length. The result is an increase in time for the ray to cross the slab. The 


opposite occurs when py < 0, which means that the ray crosses the slab in less time. That 


is, for Oo < 45°, changes in the ray path length dominate the scattering of the rays. 
A random walk model was tested. First, the sound speed fluctuations pn are 
assumed random, with (in) = 0 and (pnp) = (u? )Sam. Next, Colosi and Flatte (2007) 


derived the following: 











ul cos 8 
dx= > dx, = meses | 3.4 
sy ae J sin 6, my ee) 
coy cos26, L & 
dt=>_ dt, =t 0 — Ee , 3.5 
» ” sin’ , a sin’, Cy “a ‘ a 


where dx and dt are the summations of 5x; and dt; , respectively. Colosi and Flatte (2007) 


derived the results for the moments of dx and dt. 


(ax?) -( os ) N (42) _ 08 Ao LN(,?) (3.6) 








sin” Oo 


10 





(a) =v" 278) N (yu?) = 908 200 Lng fy) (3.7) 


sin? Oo sin® Ao co 


The growth of (dx?) and (dt*) is diffusive because they grow linearly with the number of 
scattering events (N). 


The approximations for 6x; and Oot; were numerically tested with random 
realizations and sound speedsu. The following parameter values were chosen: sound 


speed Co = 1500 ms!, depth of the layer L = 100, number of layers N = 100 and 


(8c7)*= 0.15 ms’. The variance of dx and dt were computed for an initial range of 


angles between 2 and 88 degrees. The approximations proved to be excellent, as 


demonstrated in Figure 6. 


A. DECOMPOSITION ALONG AND ACROSS THE FRONT 


The scattered rays need to be compared to the unperturbed wavefront by 
correcting the scattered ray to the same time as the unperturbed wavefront. With the 
assumption that the scattered ray propagates at the angle #0, Colosi and Flatte (2007) 


arrived at the following geometrical analysis (see Figure 7): 
Aw=dx sino, and Ar=dx cos @o-dtco (3.8) 


where Aw signifies the deviation along the wavefront and Ar is the deviation in the ray 
direction, which is normal to the wavefront. The second term in Ar represents the time 
correction back to the time of the unperturbed front. Through Equation 3.8, a critical 
understanding of the relative sizes of the two wavefront deviations is achieved. As a 


result of the fact that dx and dt have the same sign for 0, < 45°, the two contributions to 


Ar nearly cancel and the rays fall back along the unperturbed front. As a result of this 
near cancellation, Ar is reduced by a factor of tan@ relative toAw. This is a significant 
result astan@ can be quite small for launch angles near zero degrees. Glancing 
incidences created a larger distance and thus longer time for crossing the slab, 
accentuating the effect. 
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The second moments were then obtained using the random walk model described 


earlier. 





2 NI2 2 
(a2, \=NL? (WPS, (a2 eae 3.9) 


When @o approaches 90 degrees, or near normal incidence, cos” @0 approaches zero, 
thus (Aw?) becoming very small; however, sin? @0 approaches one, thus (Ar? } 
approaches a constant. As@o approaches zero or glancing incidence, (Aw?) grows more 


rapidly than (Ar”) because sin* 90 grows more rapidly thansin’ 00. At @0=45 degrees, 


the curves intersect (see Figure 8). The figure also shows that the Monte-Carlo 
simulation does very well in approximating (Ax”) and (Aw”) and is very close to 


theoretical estimates. 
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IV. 2-D CONTINUUM MODELS: NO WAVEGUIDE 


The discrete 1-D model used in the previous chapter did well to describe the 
geometry of the problem in question, but, nonetheless, brought forth serious 


disadvantages. In the 1-D model, the x-component of the ray slowness (cos Oc) was a 


conserved quantity, thus it can strongly affect the fluctuations in the system. From 
Beron-Vera et al. (2003) and Beron-Vera and Brown (2004), in ray systems in which the 
sound speed is a function of only one variable. The wave number component 
perpendicular to the sound speed gradient is constant. When the sound speed varies in 
multiple directions the wave number components are not conserved and thus, ray chaos 
can occur (i.e., exponential sensitivity to initial conditions). Further, linearization of the 
equations was prohibited in the 1-D model at an incident angle of zero as total reflection 


issues occurred. 


A. DEVIATIONS IN THE WAVEFRONT DIRECTION 


In the case of no sound channel, the geometry represents an unperturbed ray as a 
straight line from the source. The direction x is along the unperturbed ray starting at the 
point (0,0). If there are no fluctuations in sound-speed, the ray arrives at position z = 0 
and x = R after time T, as opposed to Z and R+X after time T with fluctuations ( see 
Figure 9). The deviation of Z and X is along the wavefront and perpendicular to it, 
respectively. Here, Flatte and Colosi (2007) considered an unperturbed ray angle at zero, 
thus @ will be used to symbolize the angle at position X caused by the perturbations. 
The perturbation p is responsible for the deviation in the ray angle and can be estimated 


from the ray equations. In this way, Colosi and Flatte (2007) calculated the deviation in 











angle after traveling a distance dx asd@(1—-dx, where pw is the random vertical 





gradient of fractional sound speed. Their equation for the ray angle at Range R is then: 


@(R)= |, d0=-[dxy'(x) (4.4) 


The integration path of Equation 4.4 is along the path of the unperturbed ray, 


alongside the x-axis. Thus, after a range R, the variance in the ray angle is: 
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2\_ fk R R 2\Lu _ 7 »\ LaR 
(8°) =J, dxf, dxo(u (xr) (x2) J dx (Wh?) = (Wa (4.5) 
In evaluating this integral, Colosi and Flatte (2007) made the assumption that the 
correlation length of the sound speed gradient is short in comparison to the range R in 
equation 4.5. The correlation lengths of the internal waves along and traverse to the ray 
are defined as Ly and L, respectively. Equation 4.5 shows that the rms angle with respect 


to the horizontal grows as the square root of range. 


Now a derivation for the result for the along wavefront scatter is made. Using the 
supposition that a ray propagates to some position x; and undergoes a small angular 


deflection 0 (x:) , the scattered ray will then propagate to a range R following a straight 
line and to a depth deflection of dz (R,x1) . Simple geometry (Figure 10) shows that 


25) , or = (x1) =R-x1 (4.6) 








tan O(x1)0 A(x1)= 











Next Flatte and Colosi (2007) took the sum of all deflections so, the total z-deviation at 


range R is 
(R)= * S(xi)a0- -[-(R=xi)dx. (4.7) 


The variance of Z is then 

afk R X1 X2 ' : 
(z°(R)\ =R’[ dks. an(-%!) [1-2 )(u (x1) (x2). (4.8) 
The two integration variables are within Ly of each other due to a restriction imposed by 


the correlation function of. Thus, Colosi and Flatte (2007) reduced the integral to the 


following result: 


























(z-(R)) R*[ dx (1.2) ((w)" us Ri (et) La (4.9) 


ce) Ge 


As is evident, rms deviation along the wavefront grows as the 3/2 power of the range. 
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Numerical evaluations can now be performed for @ms and Zms through the 


estimation of the internal wave factors of Equations 4.5 and 4.9. The following rms 




















values, Orns (| 6° and Zms {| 60 kmare a result of propagation to a range of 1000 km with 


sound speed fluctuation values characteristic of internal waves: Ly = 10 km, L, = 0.1 km, 


and (u) =10°. 


B. DEVIATIONS IN THE RAY DIRECTION 


The perturbed time 6T can be calculated along the ray. X is given bycodT, 
where ST is the travel-time deviation for arrival at R calculated along the unperturbed 


ray. Colosi and Flatte (2007) arrived at the following: 














X=coT(R)-R= teat 8 [, u(xz)dx (4.10) 


The variance of X is then: 


(X?)=(W?)RLa (4.11) 


1 
Comparing rms values, deviation along the wavefront (Ze)? has an rms value of 


1 
60 km in contrast to an rms value of only (x?\2 = 10 m of the deviation along the ray. 


1/2 
ie 


thousand times greater along the front than across it. 











This is a surprising result. It says that 1000, stating that deviation is one- 





c. EXPRESSIONS FOR SOUND PROPAGATION THROUGH INTERNAL 
WAVES 


The wavefront wander 1’ (Flatté et al., 1979), is directly related to the deviation of 


the wavefront along the ray (x?), As such, Colosi and Flatte (2007) derived the 


following: 
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R 
(X?) =e? =| (uw? (z)) Lo(@=0,2) dx (4.12) 
where L, is the correlation length of the internal waves in the x-direction (along the ray) 


and@. is the ray angle with respect to the horizontal. Also, Lp is dependent upon @. 


because internal waves are anisotropic. The deviation in ray angle@ is different in that it 
is not dependent upon i, but instead upon the correlation of the vertical gradient of u. 


Therefore, 
(0) =[ (uw? (z))L, (00.2) dx (4.13) 


where Ie represents the correlation length of the vertical gradient of the internal waves in 


the x-direction. Appendix A gives expressions for these correlation lengths. Colosi and 


Flatte (2007) also show 
(ZR. (1-2) (u’ (z))L, (@=0,z) dx (4.14) 


D. NUMERICAL RESULTS 


Figure 11 shows the statistics of the rms ray angle and of the rms along the 


1 1 
wavefront ‘Ze and across wavefront (x?) deviations. With regards to the range 


scaling, the rms along the front, plotted in the upper portion of Figure 10, behaves as 


1 


expected, rapidly growing as R*, Similarly, the rms ray angle (0° )? displayed in the 


1/2 


lower portion of the plot grows as expected by R’~. However, the across the wavefront 


statistics (x?) show some deviation from the prediction of Equation 3.8. Despite the 


subtraction of the parabolic curve from the Monte-Carlo simulation, at ranges beyond 


3/2 


300 km, the growth rate of (x?) ceases to increase by R'” and jumps to a rate of R*” for 


1 
reasons not known (future research). Consequently, a new approach for (x?) iS 


required for R>300 km. 
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E. OCEAN APPLICATIONS: HORIZONTAL PLANE SCATTERING 

An important result is obtained when an examination of the horizontal spreading 
of acoustic rays in the horizontal (lat, long) plane of the previous section is performed. 
The effect is that this out of plane scattering may now play an important role should the 
wavefront become obstructed by an island or seamount, which would lead to the healing 


of the wavefront. 


Colosi and Flatte (2007) modified Equation 4.9 to yield the horizontal deviation 


(Y°) for an isotropic sound speed fluctuation field: 





og: = (4.15) 


An order of magnitude can now be obtained using parameters (same as the 


previous section) typical of ocean internal waves and a propagation range of 1000 km. 














1 
The rms value of the out of plane deviation (ye) 600 m, which shows the effects of 


internal waves to appear quite small. However, for a random field of mesoscale eddies 


























with (uw?) 10° and Lu 50km, a much larger rms deviation of 10 km is obtained. 


Diffraction grows as R'” as opposed to scattering healing which grows as R°”. That is, 


diffraction dominates in shorter ranges and scattering healing dominates in larger ranges. 


F. WAVEFRONT HEALING IN THE VERTICAL 


The following images show the progression of the healing of a wavefront in the 


1 
2 


vertical past an underwater obstruction with respect to (2?) . A sound source emitting 


rays is considered. Some of the rays will hit the obstruction and others will go over it. 
As a result, a void or shadow zone will occur with respect to the vertical. As such, the 
along front scattering in question is demonstrated to reveal faster infilling of the shadow 
zones than predicted by regular diffraction (not shown). At only 100 km past the 


obstruction, the wavefront completely fills in the shadow zone. 
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V. INCORPORATION OF THE WAVEGUIDE 


The results of the no waveguide model were very encouraging, but now we want 
to treat the important effects of the waveguide. Therefore, it was necessary to modify the 
theory to incorporate the ocean waveguide and depth dependent ocean internal waves. 
To better understand this new case a geometrical approach is taken. In the no-waveguide 
case, the ray propagates to some point x’, scatters and then propagates to the final range 


R. From Snell’s law, that ray will shift some vertical differential distance d@ after 
d ; 
traveling some horizontal distance R-x’. As a result, e equals R-x’ (see Figure 10). 


The geometry in the waveguide case is similar to the no-waveguide case except that the 
ray is now following a curved ray path and not a straight line (Figure 13). Thus, the 


derivations for the variance in X, Z, and@ remain relatively the same with the exception 


that is now dependent upon the ray tube function€. Colosi and Flatte (2007) show 





that 

az». G2(x') 

zy 1 
G6) = c(R) re 


The following statistical equations from Colosi and Flatte (2007) are for the waveguide 


Case: 


(x*)=[/ ds, (w?), (@)=[, ds (uw?) L,, (Z*)=f) ds (u?) L, (5.2) 


A. NUMERICAL RESULTS 


In comparison to the results of the previous chapter, the rms ray angle 


1 
statistics (0° )’ continues to show the expected growth of R'” with some deviation at 


larger ranges followed by a tapering off at a range of about 1000 km. The across the 


1/2 


1 
wavefront statistics (X*)’ now grow as R”™ throughout the entire range interval as 
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opposed to the R*” seen in the no-waveguide case after 300 km range. Theory fits well 


for (Z? : up to 600 km, but shows deviations at greater range. For large ranges, there 


maybe a problem with the calculation of the ray tube function € past ranges of 600 km. 


Whereas, the Monte-Carlo simulation tapers off at larger ranges and grows nicely at R*?, 


the theory shows a dramatic deviation from the model statistics (Figure 14). 
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VI. CONCLUSIONS 


The ray-based physical models did very well in elucidating the physical 
mechanism for the ocean scattering processes that tend to redistribute acoustic energy 
dominantly along the wavefront. The models in question do a relatively good job 
predicting the wavefront distortions caused by the ocean internal waves with respect to 
Monte-Carlo simulations. For the simple 1-D model, dx and dt have the same sign, thus 
ray direction variance is diminished. Further, wavefront scattering is dominated by path 
length changes of the ray. In the continuum models, ray direction deviations depend on 
u and grow diffusively as R'”. Also, wavefront direction deviations depend onw and 


there is a moment arm effect that makes (2?) grow like R*?. With regards to the 


waveguide case model, the problem with the theoretical prediction of the Z statistics 
showing the dramatic deviation with the Monte-Carlo simulation must still be worked 
out. Nonetheless, the statistics provided for seven degree angle rays propagating through 
Garrett-Munk internal waves in a Munk sound speed profile remain very compelling, and 


that the distribution of acoustic energy along the wavefront is dominant. 
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VII. FUTURE APPLICATIONS 


This paper looked only at one ray with a launch angle of seven degrees using a 
ray theory approach. This method has some obvious inherent problems. First, a source 
will emit more than one ray; and, the present ray theory does not give any information 
concerning the amplitude of the acoustical field. To apply the results presented herein to 
an operational environment, the energy distribution of all the rays must be inferred by 
observing where the ray endpoints end up at along the wavefront. Further, the 
calculations should also be compared with those implementing parabolic equations and 
normal mode theory to verify results and the effects of diffraction. This new approach 
will help to understand the practical implications of along the wavefront scattering. One 
application is the rapid infilling of shadow zones behind underwater obstructions, the 
significant transfer of acoustical energy from high to low angles, and the dramatic 
vertical extension of acoustical caustic zones that abut shadow zones that were reported 
by the Acoustic Thermometry of Ocean Climate group at Navy SOSUS stations in 1999 
will be more completely elucidated. 


For wavefront healing the important result is the comparison of diffractive healing 


which scales as R!” 


and scattering which scales as R®*. At short ranges, diffraction 
dominates, but at longer ranges scattering will dominate. Where this transition occurs is 


unknown. 


Further, the theory needs to be compared to experimental data. Because a beam 
ensonifies only a small group of launch angles, an experiment involving the transmission 
of a narrow angle beam using an array of sources can meet this need. Littoral Warfare 
Advanced Development (LWAD) would be an appropriate platform to conduct this 


experiment. 
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APPENDIX A. CORRELATION LENGTHS: THE GARRETT- 
MUNK SPECTRUM 


The GM spectrum is defined by two key variables: vertical mode number j and 
frequencyo. For every vertical mode, there is some depth they are dependent upon. 


These depths are in turn locally represented by a harmonic function with a wavenumber 


given by: 
k.=ajN(z) (A.1) 
1 pD 
=— N ; 
a =a dzN(z) (A.2) 


where D is the ocean depth and N(z) is the buoyancy frequency profile. The following 


equations are expressions for the correlation length of u in the direction@, for the GM 








model: 
N 
Lo(0.2)=(7')—— °( ene) (A3) 
T Wi Wi 
= a a 
“1 ~ yt ; max 
(i [M2 oes mE) (AA) 


(A.5) 


where j,,,, 1S a cutoff mode number andqi is the inertial frequency. (u’) and N, the 


buoyancy frequency, are both functions of depth and j, is usually equal to three. 


The correlation functions for the vertical derivative of pw in the 


direction @ are obtained by giving 


L, (,z) =(aN) ( j4tn[ 8" (A.6) 


2 


1 


= A.7 
Fae a 


Imax 
(j)=M Sj 


The cutoff mode number j,,,, 18 of importance for Ly given the divergence of ( ij) 
fOr Jmax — ©- Through observation, j,,,, 18 established when the vertical wavenumber of 


the internal waves equals 0.1 cpm (Duda et al 1992). 
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APPENDIX B- FIGURES 
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Figure 1. Wavefront data collected from Navy SOSUS receivers. 
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A ray computation showing an acoustic wavefront at roughly 1,000-km range, 


resulting from a point source on the sound channel axis. The calculation was 
carried out using a Munk canonical sound-speed profile without internal-wave 


Figure 2. 


sound-speed perturbations. 
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Figure 3. Range evolution of the ray endpoints for 12,000 realizations of ray 


propagation through internal wave induced sound speed perturbations. For all 
realizations the initial angle is zero degrees and the initial depth is 1,000 m. The 
solid vertical lines show the location of the unperturbed wavefront, assumed a 
plane wave. The dash red curve shows a parabolic fit to the ray endpoint 
distributions, and the green curve shows the location of the unperturbed point 
source front. 
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Figure 4. Range evolution of the ray endpoints for 12,000 realizations of ray 
propagation through internal wave induced sound speed perturbations, with the 
Munk canonical profile. For all realizations the initial angle is 7 degrees and the 
initial depth is 1000-m. The red curve shows the unperturbed wavefront from a 
point source, and the green dot shows the location of the unperturbed 7 degree 
ray. 
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Unperturbed Ray 








Figure 5. Schematic showing the geometry of the 1-dimensional discrete layer problem. 
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<dyx2>' ys, 6, 
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<dx*> theo 


<dt2> calc 


<dt2> theo 





Figure 6. Comparison of statistics for dx and dt as a function of the initial angle @o, 
using the random walk model. The blue curves represent the Monte Carlo 
simulations and the red curves are the theoretical estimations. 
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Figure 7. Schematic of the geometry of the projection of the scattered ray onto the 


unperturbed wavefront. 
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Figure 8. Computations of (Aw?) and (Ar”) as a function of initial o, from direct 


Monte-Carlo numerical simulation and from theoretical estimates using Eq. 3.9. 
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Figure 9. Geometry of ray scattering in the absence of a waveguide. 
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Figure 10. Geometrical description of the application of Snell’s Law to the no-waveguide 


Case. 
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Figure 11. Computations of rms ray deflections for zero angle rays propagating through 
GM internal waves (solid) and theoretical expectations (dash). The rms angle is 
diplayed in the lower portion of the plot, showing the expected growth by R? 
The rms along the wavefront deflection is diplayed in the upper portion of the 

plot, also showing the expected Re growth. The ray direction deflection is 

displayed in the middle section of the plot where the numerical result has been 
corrected for the parabolic spreading of the ray endpoints. The R'? scaling breaks 

down after 300-km in range. 
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Figure 12. Progression of along the wavefront scattering past an underwater obstruction 


demonstrating a complete infilling of the shadow zone at a range of 100 km. 
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Zz scattering event 


Figure 13. Geometrical description of the application of Snell’s Law to the waveguide 
case. 
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Figure 14. Computations of rms ray deflections for 7 degree angle rays 
propagatingthrough GM internal waves in a Munk sound speed profile. Monte- 
Carlo simulation results are shown with solid lines while theoretical results are 

dashed. The rms ray angle is in the lower portion of the plot and the expected 
display of R’” is evident. The rms along wavefront deflection is displayed in the 
upper portion and the expected growth of R3/2 is evident. The rms across the 


wavefront is in the middle of the plot and shows the expected growth of RY, 
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